Chirality-dependent phonon-limited resistivity in multiple layers of graphene 
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We develop a theory for the temperature and density dependence of phonon-limited resistivity 
p(T) in bilayer and multilayer graphene, and compare with the corresponding monolayer result. For 
the unscreened case, we find p « CT with C oc «p 2 in the high-temperature limit, and p w AT 4, 
with A oc Vp 2 fcp 3 in the low-temperature Bloch-Griineisen limit, where vf and kp are Fermi velocity 
and Fermi wavevector, respectively. If screening effects are taken into account, p ~ CT in the high- 
temperature limit with a renormalized C which is a function of the screening length, and p ~ AT 6 
in the low-temperature limit with A oc fep 5 but independent of vf- These relations hold in general 
with vf and a chiral factor in C determined by the specific chiral band structure for a given density. 
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The hallmark, indeed the universally used definition, of 
a metal is its phonon-scattering induced resistivity with 
increasing temperature. Understanding, and if possible, 
controlling the electron-phonon interaction is thus one 
of the most important fundamental physical problems in 
any new electronic material. Since electron-phonon cou- 
pling typically controls the room-temperature resistiv- 
ity of all metals (and doped semiconductors), the study 
of phonon-limited resistivity is also crucial for techno- 
logical applications. In this Letter we study theoreti- 
cally the electron-phonon interaction induced resistivity 
in graphene of arbitrary layer thickness (including mono- 
layer, bilayer, and multilayer graphene), finding a num- 
ber of new experimentally testable results of considerable 
importance. In particular, we predict that the exponent a 
in phonon-induced graphene resistivity, p ~ T a , depends 
on the characteristic of screening and at low enough tem- 
peratures the resistivity eventually goes as p ~ T 6 for 
screened phonon scattering, and that chirality has a sub- 
tle quantitative effect on the high-temperature linear-in- 
T resistivity. 

When the phonon energy (tvjJq) is much lower than 
the Fermi energy (ef), he. huj q -C ef, the scattering of 
electrons from acoustic phonons can be divided by two 
regimes: T < Tbg an d T > Tbg- The characteristic 
temperature Tbg is known as the Bloch-Griineisen (BG) 
temperature and is defined as k^T^G = 2hv p i l k F , where 
kv and i> p h are the Fermi wavevector and the sound veloc- 
ity, respectively [l[ . For T > Tbg > the number of phonons 
increases linearly with temperature and so does the resis- 
tivity limited by electron-phonon scattering, p oc T. This 
behavior in graphene has been observed experimentally 
0-13] • In theoretical work [j| it is found that the temper- 
ature dependent resistivity of monolayer graphene in BG 
regime is given by p oc T 4 in the absence of screening. In 
a recent careful measurement of the temperature depen- 
dent resistivity of a high density graphene @ , a smooth 
transition of the resistivity from a linear T dependence 
to a T 4 dependence is observed as the temperature de- 
creases below Tbg- The measured resistivity in Ref. @ 
is the first explicit observation of BG behavior in two- 



dimensional (2D) systems. Graphene may be an ideal 
system to observe the BG behavior because at relatively 
high densities (n > 10 13 cm" 2 ), T BG > 200 K and the 
other extrinsic scatterings are severely suppressed [l| ■ In 
our current work, we study a general multilayer graphene 
system including screening effects of the electron-phonon 
interaction, and demonstrate that the power-law depen- 
dence could vary smoothly from a T 4 to T 6 power law 
depending on the screening strength. 

Motivated by the recent measurement of BG behav- 
ior of phonon limited resistivity 0] we investigate the 
intrinsic transport properties of bilayer and multilayer 
graphene as limited by phonon scattering using Boltz- 
mann transport theory in the case of &bT <C £f- he. 
high-density systems. In this paper, we consider only 
the longitudinal acoustic phonons because other phonon 
modes are negligible in the temperature range of our in- 
terest [f|. 

From the Boltzmann transport theory, the energy- 
averaged relaxation time in k^T -C e F limit is given by 
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where 2/(ef) is the density of states per spin and val- 
ley at the Fermi energy e f , |C(fc F )| 2 = is the 
squared matrix element for acoustic phonon scattering 
at the Fermi wavevector fcp, D is the acoustic phonon 
deformation potential, p m is the graphene mass density 
and v p h is the phonon velocity. The integration factor / 
in Eq. (flj is given by 
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where q = 2/cf sin(<^>/2) is the magnitude of an acous- 
tic phonon wavevector, F(q) is the chiral factor defined 
by the square of the wavefunction projection between in- 



coming and scattered states, N q 
phonon occupation number, uj q 



(e^« - I)- 1 is the 
i» p h<z is the acoustic 



2 



phonon angular frequency, e(q) is a dielectric function 
and (3 = l/(k B T). 

Then, the conductivity in 2D system is given by 
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and corresponding resistivity is p — a 1 , where v-p is the 
Fermi velocity defined by wp — j^\ £ _ Ef , 9s and g v are 
spin and valley degeneracy factors, respectively. When 
several bands cross the Fermi energy, we add conductivity 
contributions from each crossing band (with the same 
screening wavevector determined from the total density 
of states). 

The dielectric function e(q) takes into account the 
screening effect at wavevector q. Within random phase 
approximation, e(q) — 1 + q s (q)/q, where q s {q) is the 
screening wavelength [8]. In our temperature range, 
we can approximate q s (q) ~ <7tf where ^tf is the 2D 
Thomas-Fermi screening wavevector defined by ^tf = 
g s g v cigr(v / vp)k-F where a gr = e 2 /(ehv), e is the effec- 
tive dielectric constant and v is the monolayer in-plane 
velocity Q- In 2D, the strength of screening is de- 
termined by the parameter qo = ^tf/^f oc £% r Q> 
thus unscreened (strong screening) limit corresponds to 
qo,a gI -> (q ,a gI > I). 

By setting x = q/(2k F ) = sin(</>/2), Eq. © can be 
reduced to 
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where xtf = <Ztf/(2£;f) and Zbg = Tqg/T. 

In the high-temperature limit (T > T B g), Eq. (0} be- 
comes I s» Zoo/zbg where 
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thus the resistivity becomes p ~ CT where 
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Note that in the high-temperature limit, all phonons are 
thermally excited giving the linear T dependence of the 
resistivity. Eq. ([5]) contains the chiral factor F, thus C 
depends on the chiral properties of wavefunctions. The 
screening effect enters only in the integration factor of 
Zoo, thus does not change the temperature dependence 
qualitatively. Note that in the unscreened case (a gr = 0), 
C oc Up 2 , while in the strong screening limit (a gr 3> 1), 
Zoo oc Vp, thus C oc v F , independent of up and density. 

Next, we consider the low temperature limit (T <C 
Tbg)- In the unscreened case, after setting y = z B gx, 
Eq. (@} becomes 
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FIG. 1. Acoustic phonon-limited resistivity of monolayer 
graphene as a function of temperature for (a) o% r = with 
several densities, (b) a g 



n = 10 13 cm" 2 
derivatives of (c" 



1 with several densities, (c) 
with several a gr , and (d) the logarithmic 



where £(s) is the Riemann-zeta function and F(0) = 1 
was used. Thus, the resistivity becomes p « AT 4 " where 
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Note that A is independent of the chiral factor F due 
to the small angle electron-phonon scattering. The T 4 
dependence of the resistivity arises from phase space 
limitations for electron-phonon scattering as phonons of 
wavevector 2k F are no longer thermally excited. The 
transition from T and T 4 dependence occurs around 
T a Tbg- 

In the screened case, Eq. iQ becomes 
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thus the resistivity in the low-temperature limit becomes 
p « AT 6 where 



A = 
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Note that compared to the unscreened case in Eq. (jHJ , A 
depends on k F but not on v F , and is thus independent of 
a specific band structure for a given k F or density. 

The theory developed so far is valid for a general 2D 
electron system in the case of k^T <C £ F , thus if we know 
the Fermi velocity v F and the chiral factor F(q), analytic 
expressions of the coefficients A and C can be obtained 
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from Eqs. (|6]), ([8]) and ((10)) . Note that for monolayer 
graphene, wp = v and F(q) = (1 + cos</>)/2 with q = 
2fcp sin((^)/2), whereas for bilayer graphene, 



hvkv, 



Vp 



v /(ii/2) 2 + (hvk F ) 2 
and the chiral factor for low energy band is Q 



(11) 



F(q) = r?+(l + r?) cos 0] 2 , 



where 77 = + n/no 

fiufco = t±/2. 

Thus, for monolayer graphene, we get 
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A L G(a g r 7^ 0) = A MLG (a gr 7^ 0), 

where MLG and BLG stand for monolayer graphene and 
bilayer graphene, respectively. 

Note that in the low-density limit (n <C no or 77 w 1), 
it can be shown that the 4-band bilayer results obtained 
here approach the 2-band low-energy bilayer results. In 
the high-density limit (n > no or ij w 0), 4-band bilayer 
A BLG approaches the monolayer A MLG , while correspond- 
ing C B lg does not approach the monolayer C ML g ■ 

The reason for the discrepancy in C between mono- 
layer and 4-band bilayer results in the high-density limit 
is that in bilayer graphene, the chiral factor with the in- 
terlayer hopping tj_ = cannot be obtained from the 
limit t±_ — > because wavefunctions involved in the two 
cases are different, though the energy spectra become 
similar. Note that the coefficient A is independent of 
the chiral factor due to the small angle scattering, thus 
monolayer and bilayer results converge in high enough 
density. 

For numerical calculations, we use p m — 7.6 x 1CP 8 
g/cm 2 , u p h = 2.6 x 10 6 cm/s and D = 25 eV for phonons, 
and the nearest-intralayer hopping t = 3 eV and nearest- 
interlayer hopping t± = 0.3 eV. For simplicity, other 
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FIG. 2. Acoustic phonon-limited resistivity of bilayer 
graphene as a function of temperature for (a) a gr = with 
several densities, (b) a gr = 1 with several densities, (c) 
n — 10 13 cm -2 with several a gr , and (d) the logarithmic 
derivatives of (c). 



remote hopping terms are neglected keeping rotational 
symmetry in the energy spectrum. 

Figures [lja) and (b) show the longitudinal acoustic 
phonon-limited resistivity of monolayer graphene as a 
function of temperature for unscreened case (a gr = 0) 
and screened case with Q?g r — 1, respectively, for sev- 
eral densities. For the unscreened case, the resistivity 
increases as p ~ T 4 at low temperatures and p ~ T at 
high temperatures, while for the screened case, p ~ T e at 
low temperatures and p ~ T at high temperatures. The 
transition occurs around Ibg, which is given by 70.4, 
222.6 and 583.6 K for n = 10 12 , 10 13 and 10 14 cm" 2 , re- 
spectively. FiguresQJc) and (d) show the temperature de- 
pendence of resistivity at n = 10 13 cm -2 for different a gr 
values and their logarithmic derivatives, respectively. As 
a gr increases from 0, the low temperature power changes 
from 4 and approaches to 6, as shown in Fig. QJd). Bi- 
layer graphene results show very similar behavior, as il- 
lustrated in Fig. O except that the density dependence 
of the coefficients A and C are qualitatively different from 
the monolayer expected from Eq. (fT4]) . 

Let us now consider the density dependence of the coef- 
ficients A and C in a general multilayer graphene system. 
At low energies or equivalently at low densities, arbitrar- 
ily stacked multilayer graphene is described by a super- 
position of pseudospin doublets 10], as summarized in 
Tab. |TJ Thus each pseudospin doublet contributes to the 
density dependence of the coefficients at low densities. At 
high energies or equivalently at high densities, interlayer 
coupling becomes negligible and energy band structure 
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TABLE I. Chirality decomposition f or p eriodic AB and ABC 
stacking up to N = 4 layer stacks [10(. Here we have arbi- 
trarily labeled the first two layers starting from the bottom 
as A and B. 
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FIG. 3. Density dependence of (a) A for a gr = 0, (b) A for 
a gr = 1, (c) C for Q gr = 0, and (d) C for a gr = 1 with 
different stacking sequences. 



looks like that of monolayer graphene, thus the density 
dependence will follow that of monolayer graphene. 

First, consider the density dependence of A at low den- 
sities (but still assuming in k-gT <C e-p limit). Note that 
for J-chiral system, Sk oc k J , thus vf oc /c f _1 . From 
Eq. © and Eq. (fTU|) . 



A(a gT = 0) oc w F 2 fc F 3 oc n 
A(a gv ± 0) oc fcf 5 oc n~ 5/2 . 



-(J+l/2) 
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Thus, at low densities, for the unscreened case (a gr = 
0) , A has different density dependence depending on the 
chirality while for the screened case (cv gr ^ 0), A has the 
same density dependence, irrespective of the chirality. 

Figure [3](a) and (b) show the density dependence of 
the coefficient A for the unscreened case and screened 
case with a gr — 1, respectively, for periodic stackings up 
to N = 4 layer stacks. The numerical results agree with 
the preceding analysis. 

Next, consider the density dependence of C. From 



C(a gr = 0) oc «p 2 oc n- {J - l \ 
C(a gI > 1) oc v F oc n°. 



(16) 



Thus, at low densities, for the unscreened case (a gr = 0), 
C has different density dependence depending on the chi- 
rality while in the strong screening limit (a gr 3> 1), C is 
independent of density, irrespective of the chirality. Note 
that for monolayer graphene, C remains constant with 
density for both unscreened and screened cases because 
v-p = v is constant in Eq. ([5]). 

Figure [3](c) and (d) show the density dependence of 
the coefficient C for the unscreened case and screened 
case with ev gr = 1, respectively, for periodic stackings 
up to N = 4 layer stacks. The numerical results agree 
with the preceding analysis, and at high densities density 
dependence of C becomes weaker following the density 
dependence of monolayer, even though a gr is not in the 
strong screening limit. Note kink structures appear at a 
density when multiple bands begin to contribute in a mul- 
tilayer stack. Also note that for ABA and ABAB stack- 
ings, multiple bands contribute to the resistivity even at 
low densities, thus they show relatively different behav- 
ior compared with other stackings due to the combined 
effects of the chirality and multiband screening. 

Even though the phonon parameters of graphene are 
well defined, there is uncertainty about the value of the 
deformation potential D 0, 0, [ll|. The results of this 
paper can be applied to extract the proper value of the 
graphene deformation potential. To get the correct value, 
however, it is crucial to know the significance of screen- 
ing in electron-phonon scattering. Once screening is in- 
cluded, the larger value for D is required to match the 
result compared with that without screening. Since two 
approaches give rise to the same linear-in-T resistivity 
at high temperatures, it is hard to get the correct value 
of deformation potential by considering only the high- 
temperature result. 

At low temperatures (BG regime), however, the tem- 
perature dependence of the resistivity strongly depends 
on screening, i.e. ev gr , and the inclusion of screening 
modifies the behavior of p(T) from T 4 to T 6 . In this 
sense, the recent experimental observation of resistivity 
at BG regime, p oc T 4 around 10 K 0] indicates su- 
perficially that screening in the effective electron-phonon 
interaction is not strong enough to show p oc T 6 depen- 
dence. It is possible, however, that at lower temperatures 
the power-law behavior eventually becomes p ~ T 6 , as 
demonstrated in Figs, [ljd) or HJd), since the screened 
coupling should be the appropriate interaction in an elec- 
tron liquid. We suggest that the observed BG temper- 
ature exponent of 4 in Ref. 6; is simply an effective ex- 
ponent which will increase to an exponent of 6 in the 
asymptotic T <C ?bg temperature regime. For the ex- 
perimental observation of this prediction, it is important 
to work in a regime where the graphene temperature de- 
pendence from non-acoustic-phonon mechanisms [12| is 
quantitatively unimportant or can be subtracted out in 
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an unambiguous manner, which may necessitate working Authors thank D. K. Efetov and P. Kim for helpful 
at high carrier density and on special substrates, as has discussion. The work is supported by the NRI-SWAN 
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